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ABSTRACT 

We present 7-ray observations with the Large Area Telescope on board the Fermi 
Gamma-Ray Telescope of the nearby radio galaxy Centaurus A. The previous EGRET 
detection is confirmed, and the localization is improved using data from the first 10 
months of Fermi science operation. In previous work, we presented the detection of the 
lobes by the LAT; in this work, we concentrate on the 7-ray core of Cen A. Flux levels 
as seen by the LAT are not significantly different from that found by EGRET, nor is 
the extremely soft LAT spectrum (F = 2.67 ± QAOgtat i 0.08^^5 where the photon flux 
is ^> cx E^^). The LAT core spectrum, extrapolated to higher energies, is marginally 
consistent with the non-simultaneous HESS spectrum of the source. The LAT obser- 
vations are complemented by simultaneous observations from Suzaku, the Swift Burst 
Alert Telescope and X-ray Telescope, and radio observations with the Tracking Active 
Galactic Nuclei with Austral Milliarcsecond Interferometry (TANAMI) program, along 
with a variety of non-simultaneous archival data from a variety of instruments and wave- 
lengths to produce a spectral energy distribution (SED). We fit this broadband data set 
with a single-zone synchrotron/synchrotron self-Compton model, which describes the 
radio through GeV emission well, but fails to account for the non-simultaneous higher 
energy TeV emission observed by HESS from 2004-2008. The fit requires a low Doppler 
factor, in contrast to BL Lacs which generally require larger values to fit their broad- 
band SEDs. This indicates the 7-ray emission originates from a slower region than that 
from BL Lacs, consistent with previous modeling results from Cen A. This slower region 
could be a slower moving layer around a fast spine, or a slower region farther out from 
the black hole in a decelerating flow. The fit parameters are also consistent with Cen 
A being able to accelerate ultra-high energy cosmic-rays, as hinted at by results from 
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the Auger observatory. 

Subject headings: galaxies: active — galaxies: individual (Centaurus A) — galaxies: 
jets — gamma rays: galaxies — radiation mechanisms: nonthermal 



1. Introduction 



41 Radio galaxies exhibiting jets which terminate in r adio lobes on tens of kp c to Mpc scales are 

42 classified based on their radio morphology and power bv lFanaroff RilevI (|1974| ). They are divided 

43 into Fanaroff-Riley (FR) type I and type II, where type I sources have the highest surface brightness 

44 feature at the center, while in type II sources it is farther from the core. Furthermore, the transition 

45 radio luminosity betwe en FRI and FRII increases with the optical luminosity of the host galaxy 

46 (jLedlow fc Owenlll996l ). In the AGN unification scheme, blazars are thought to be radio galaxies 

47 with the jet aligned along our line of sight, and are subdivided into flat spectrum radio quasars 

48 (FSRQs) and BL Lacertae objects based on the strength of emission lines in the ir spectrum, where 



49 FSRQs generally have strong emission lines , while BL Lacs have weak or none (jStrittmatter et al 



1972: 


Marcha et al. 


1996: 


Landt et al. 


2004) 



51 BL Lacs, while FRIIs correspond to misaligned FSRQs (e.g. JUrrv &i Padovanilll995l. and references 

52 therein), although there is evidence that this unification scheme is too simple (e.g.. lLandt &: Bignall 



20081 ). Apparent superluminal motion observed on milli-arcsecond size scales indicates that their 



54 jets must be moving at high relativistic speeds, with b ulk Lorentz factor F, ~ 10-20 for FSRQs 



55 and BL Lacs (IKellermann et al.ll2004l : iLister et al.ll2009l ). although some TeV BL Lacs have Tj ~ 3 



56 ( Finer et al. 20081 ). The existence of high energy and very high energy (VHE) 7-rays observed 

57 from these sources provides further evidence for highly relativi stic flows, as they are ne cessary to 

58 avoid 7-ray attenuation by electron-positron pair production (jPondi &: Ghisellini Indeed, 

59 this sometimes gives values of Tj greater than that found from very-long baseline interferometry 

60 (VLBI) s uperluminal observati o ns; e.g., F, > 50 is required for a recent outburst from PKS 2155- 

61 304 (e.g., iBegelman et"al]l2008l : iFinke et ahlbood ). 



Since blazars are strong sources of beamed 7-rays, it is natur al to think that radio g alaxies may 



63 be also. Several radio g alaxies were detected by EGRET: 3C 111 (^Hartman et al.ll2008l). NGC 6251 



64 (jMukherjee et ahlbood ). and Centaurus (Cen) A (jSreekumar et al.lll99d : lllartman et al.1 19991 ). The 



65 identiflcations were rather uncertain, due to the large EGRET error circles. Only two radio galaxies 

66 have been detected so far with the latest generation of TeV atmospheric Cherenkoy telescopes. 



67 M87 (Aharonian et al 



20061 : lAcciari et al.1 12OO8I : I Albert et al.ll2008l : lAcciari et al.ll2009l ) and Cen A 



Aharonian et al.ll2009l ). The Radio Galaxy 3C 66B seems to have been seen by MAGIC lAliu et al 



20091 ). although the detection i s questionable due t o its proximity to the BL Lac 3C 66A and its 



71 the detections of N GC 1275 (Per A; 



lack of detection by VERITAS (lAcciari et al.ll2009l). Th e Fermi -LAT collaboration has reported 



Abdo et al. 


2009b). M87 ( 


Abdo et al. 


2009d) 



72 (jAbdo et al.ll2009d ). Several more 7-ray detections of radio galaxies have been reported in the flrst 
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73 Fermi-L AT catalog (IFGL; Abdo et al. 2010al lbl) and a future publication will examine them in 

74 more detail (Fermi Collaboration 2010, in preparation). 

75 The Fermi Gamma Ray Space Telescope was launched on 2008 June 11 and contains the Large 

76 Area Telescope (LAT), a pair co nversion telescope w hich has a field of view of about 20% of the 

77 sky at 20 MeV to over 300 GeV (lAtwood et alJl2009l ). For the first year of operation, Fermi was 

78 operated in a sky-survey observing mode, wherein the LAT sees every point on the sky every 3 

79 hours. 

80 During the first 3-months of science operation, the Fermi-h AT confirmed ( Abdo et al. 2009al l3) 

81 the EGRET detection of Cen A. Here with additional monitoring, we present accumulated data 

82 after 10 months of operation. The new LAT observations bridge the gap between EGRET and 

83 HESS, providing a detailed look at the 7-ray spectrum essential for addressing emission models. In 

84 addition to the LAT 7-ray source in the central few kpc (hereafter the 7-ray "core"), 7-rays from the 

85 giant lobes of Cen A have also been seen with Fermi, with the origin likely to be Compton scattering 

86 of the cosmic micr owave backgrou nd ( CMB) and extraga l actic background light (EBL), confirming 

87 the predictions of ICheund (j2007l ) and iHardcastle et al.l (120091). Detailed work on separating the 

88 core and lobe emission is presented elsewhere (jAbdo et al.ll2010d . hereafter referred to as the lobe 

89 paper), although we provide a summary of LAT observations below. For the purposes of this paper, 

90 which is a study of 7-ray emission of the core, the lobes are essentially background sources. 

91 We present a summary of Cen A and observations of this object in § [2l The observations of 

92 the core of Cen A with the LAT over the first 10 months of Fermi operation are presented in § [3j 

93 We also present simultaneous Cen A core observations from Suzaku and Swift, and radio data from 

94 the TAN AMI program in § [H In § [5] we combine these with archival data and model its SED of the 

95 Cen A core. In § [6] we discuss the implications in detail, and we conclude with a brief summary (§ 

96 El). 



Centaur us A 



98 The FRI Cen A is the nearest radio found active galaxy to Earth, making it an excellent 

99 source for studying the physics of relativistic outflows and radio lobes. Indeed, it is near enough 

100 that its peculiar velocity dominates over the Hubble flow, and its redshift {z = 0.00183) cannot 

101 be used to accurately calculate its distance. iFerrarese et al.l (j2007l ) have found that the average 

102 of several distance indicators gives D = 3.7 Mpc, which we adopt. At this distance, an arcsecond 

103 corresponds to about 18 pc. Due to its proximity to Earth, it has been well studied throughout 

104 the electromagnetic spectrum, from radio to 7-rays. Recently, the Auger collaboration reported 

105 that the arrival directions of the highest energy cosmic rays (> 6 x 10^^ eV) observed by the Auger 

106 observatory are correlate d with nearby AGN, including Cen A ([Abraham et al.l 120071 . l2008l ). while 



Moskalenko et al.l (|2009l ) found that, if the giant lobes are taken into account, as many as four 



108 ultra-high energy cosmic rays (UHECRs) may be associated with this source. Although the overall 
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109 significance of this correlation i s reduced in tlie expa nded Auger data set, tlie significance remains 

110 liigh in the direction of Cen A (jAbraham et al.ll2009l ). This suggests that Cen A — and other radio 

111 galaxies — may be sources of UHECRs. 

112 Cen A has interesting radio structure on several size scales. The most prominent features 

113 are its giant radio lobes, which subtend ~ 10° on the sky, oriented prim arily in the North-S outh 

114 direction. They have been im aged at 4.8 GHz by th e Parkes telescope (jJunkes et al.l Il993l ) and 

115 studied at up t o ^60 GHz by 



116 Probe (WMAP; iHinshaw et al 



Hardcastle et al.l (|2009l ) utilizing Wilkinson Microwave Anisotropy 



2009 ) observations. The North lobe co ntains a bright region a few 



117 tens of arcminutes in size often referred to as the Northern middle lobe (jMorganti et al.lll999l ) . Mis- 

118 aligned by approxi mately 45° relative to the outer lobes are inner radio lobes on an arcminute scale 

119 ( Burns et al A strong, well-collimated jet can be seen on the arcsecond size scale in the radio, 

120 and Chandra can resolve X-ray eniission from it, which is likely caused by synchrotron emission 

121 ([Kraft et al.ll2002l : iHardcastle et al.ll2003l ). The inn ermost region of Cen A has been resolved wit h 



122 VLBI, and shown to have a size of ~ 3 x 10^^ cm (iKellermann et al 



1997 



Horuichi et al.l 120061 ) . 



123 Observations at shorter wavelengths also re veal a small core, namely VLT infrared interferometry 

124 which resolves the core size to ~ 6 x 10^^ cm ( Meisenheimer et al.ll2007l l. VLBI images reveal a weak 

125 counter jet on the milli-arcsecond scale (| Jones et al.lll996l ). Based on the motion of th e VLBI blobs. 



126 and as suming the brightness differences of the different jets are due to Doppler effects, iTingav et al 

127 (jl998l ) estimate the angle of the sub-parsec jet to our line of sight to be ~ 50 — 80°. Applying a 



28 simila r technique to the 100 pc scale jet which has a larger j et-counterjet ratio, IHardcastle et al 

29 (|2003l ) estimate a jet angle of ~ 15°. IHardcastle et al.l (120031 ) speculate that the conflicting angle 



30 estimates may be due to the assumption that the jet-counter jet brightness differences are caused 

31 by Doppler beaming rather than intrinsic differences. 

32 NGC 5128, the giant elliptical host galaxy of Cen A, contains a kiloparsec-scale dust lane. This 

33 feature appears to be an edge-on di sk obscuring the ce n tral region and nucleus, and is probably 



34 the remnant of a previous merger (jQuillen et al. 



1992 



Israeli Il998l ) . It also has a dusty torus 



35 withi i i 100 pc of the b lack hole, with a high column density {Nh > 10 cm ) (llsrael et al 



20081 : IWeifi et al.ll2008l ). X-ray spectra taken at various times over decade timescales indicate a 



37 time-var ying absorbing column , density, which could be due to variations in a warped disk viewed 

38 edge-on (jRothschild et al.ll2006l ). Estii nates for the mass of the supermassive black hole at the center 

39 of Cen A range from (0.5-1) x 10^ Mq (jsilge et ahlboOsUMarconi et al.ll2006l : lNeumayer et al.ll2007l ) 

40 based on the kinematics of stars, as well as H2 and ionized gas. 



W ith the Compton Gamma- Ray Obse rvatory, emission was detected by OSSE (iKinzer et al . 
1993) and COMPTEL dsteinle et alJll998l l at 100s of keV to MeV energies. Ixinzer et al.l ^199^ ) 



43 suggested the hard X-ray emission from Cen A detected with OSSE was the result of Compton- 
.44 scattered disk radiation by a therr nal plasma (i. e ., a ho t corona), due to a turnover in the spectrum 

45 at a few hundred keV. However, ISteinle et al.l (| 19981 ) noted that the high-energy portion of the 

46 OSSE spectra smoothly connected with the higher energy COMPTEL spectra, and the OSSE and 
COMPTEL variability seem to be correlated. They used this to argue for a nonthermal jet origin for 
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148 the X-rays. Evans et al.l (|200J) have resolved the arcsecond-scale core of Cen A with Chandra and 

149 XMM-Newton. The 2-7 keV X-ray continuum, when corrected for absorption, is consistent with 

150 what is predicted from a, correlation betwee n unresolved X-ray emission and 5 GHz core emission 

151 for jets of radio galaxies (jCanosa et al.lll999l ). They thus consider it likely that nonthermal emission 

152 from the sub-pc (sub-mas) scale jet is the origin of the continuur n X-rays from t he co re of Cen A. 

153 However, hard X-rays observed by Su zaku do not seem to fit on the lCanosa et alJ (jl999l ) correlation, 

154 possibly indicating a non-jet origin (jMarkowitz et alJ 120071 ) . The nature of the continuum X-ray 

155 emission from the core of Cen A remains an open question. 



156^ Cen A has been a target of 7-ray observations dating back to the 197 0s fe.g.. iGrindlay et al 



1975 




Hall et al. 


1976) 


Hartman et al. 




1999^. 



Hall et al.l 1 19761 ). Cen A was seen by EGRET up to GeV energies (jSreekumar et al.lll999l : 
I999I ). The 7-rays are thought to originate from a relativistic jet near the cen- 

159 tral elliptical galaxy (the radio "core" ) analogous to blazars, although it has been suggested that 

160 Compton-scattering of the C MB and the infrared-optical EBL in t he giant radio lobes could be a 

161 source of 7-r avs from Cen A (IHardcastle et al.ll2009l : ICheungj 120071 ) and other radio galaxies such 

162 as Fornax A (jGeorganopoulos et al.ll2008l ). At the highest, TeV energies, a detect ion was recently 

163 reported from Cen A by the air Cherenkov detector HESS (jAharonian et al.ll2009l ). 



3. FermirhAT Gamma-Ray Observations 
3.1. Localization 



166 The EGRET detection of Cen A (jSreekumar et al.lll999l : iHartman et al.lll999l ) was confirmed 

167 early on by the Fermi-L AT . Based on 3- months of all-sky survey data , the in itial LAT detection 

168 was reported in the LAT bright source list (BSL) paper (jAbdo et al.l l2009al ) as OFGL J1325.4- 

169 4303 with a 95% confi dence localization 95=0.304°= 18.3'. In the companion LAT Bright AGN 

170 Sample paper (LBAS: lAbdo et al.ll2009a ) to the BSL, a single power-law fit was reported, which 

171 gave F(>100 MeV) = 2.15 (± 0.45) x lO"'^ ph cm'^ s'^ with photon index, T = 2.91 ± 0.18, 

172 and a peak flux on a ~1 week timescale of (3.23 it 0.80) x 10^^ ph cm^^ s^^. Note that this only 

173 considered the 7-ray emission from Cen A as a single point source, i.e., it did not account for any 

174 lobe emission. 

175 To these initial observations, 7 additional months of all-sky survey data are added to the 

176 current analysis. Specifically, the observations span the time period from 2008 August 4 to 2009 

177 May 31, corresponding to MET (mission elapsed time) 239557420 - 265507200. Diffuse event class 

178 (CTBCLASSLEVEL=3) events were selected with a zenith angle cut of <105°, and a rocking angle 

179 cut of 39°. The former are well calibrated and have minimal background while the latter greatly 

180 reduce Earth albedo 7-rays. For the analysis, LAT Science Tool^ version v9rll was utilized with the 



^http : //f ermi . gsf c . nasa . gov/ ssc/data/analysis/scitools/ overview . html 
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181 P6_V3_DIFFUSE instrument response function (IRF). The standard LAT Galactic emission model, 

182 GLL_IEM_V02 . FIlH was used and the uniform background was represented by the isotro pic diffuse 7- 

183 ray ba ckground and the instrumental residual background 



2009a). We consider 11 point sources in the IFGL catalog (jAbdo et al.ll2010al . see also Figured]) 



isotropic ieni_v0 2 . txt , 



Abdo et al, 



185 Figure [U shows the the 0.2—30 GeV LAT image centered on Cen A, which is clearly detected. 

186 Also prominent is the Galactic emission toward the south, and several faint sources in the field. We 

187 obtained a localization of the source at Cen A with gtf indsrc, which finds point source locations 

188 based on an unbinned likelihood analysis. The resulting localization was reduced to 95= 0.087° 

189 = 5.2' (5.7 kpc), centered at RA = 201.399°, Dec = -43.033° ( J2000.0 epoch) which is offset by 

190 0.029° = 1.7' (1.9 kpc) from the VLBI radio position of Cen A (|Ma et ahlEggsl l. Figure [2] shows 

191 the localization error circle of the LAT emission overlaid on the combine d radio, optical, X-ray 



192 images. The new LA T position is consistent with that of 3EG J1324-4314 (jSreekumar et al 



Hartman et al 



199S 



19991'! ■ but bo t h are notablv offset from EGR J1328-4337, the closest EGRET source 



194 in the lCasandjian &: Grenied (j2008l ) catalog. The latter derived position shifted in such a way that 

195 Cen A was outside of the 95 localization circle, so that there was some ambiguity as to whether 

196 EGRET was actually detecting Cen A, but the new LAT position confirms the earlier 3EG result. 

197 The LAT significantly improves upon the previous EGRET 7-ray localization (95=0.53°= 32'). 



3.2. Spatial and Spectral Analysis 



199 The binned likelihood fitting was performed with the gtlike tool, first assuming Cen A is a 

200 point source, i.e., that there is no 7-ray lobe emission (model A). The field point source positions 

201 were fixed, and their spectra were assumed to be power-laws, with the photon indi ces allowed to 

202 vary. The location of Cen A was fixed at its VLBI radio position (iMa et al.lll998l ). In addition 

203 to the 11 IFGL point sources used in the lobe paper, in order to treat the lobe emission as a 

204 background source, we include two IFGL sources, IFGL J1322.0 - 4515 and IFGL J1333.4 - 4036, 

205 which are thought to be the local maxima of the lobe emission. A likelihood analysis with the energy 

206 information binned logarithmically in 20 bins in the 0.2-30 GeV band, and the 7-ray directions 

207 binned into a 14° x 14° grid with a bin size of 0.1° x 0.1°. For both the Galactic and isotropic 

208 emission models, one free parameter was introduced to adjust the normalization. Because the 

209 effective area of the LAT is rapidly changing below ~ 200 MeV, we use events with energy above 

210 this value. Above 30 GeV the significance of detection is < 3a, so we make a cut as this energy as 

211 well. 



212 As a result, the test statistic (TS; iMattox et al.l 119961) is found to be 378 for Cen A, which 

213 is smaller than the TS=628 in the IFGL catalog (jAbdo et al.ll2010al ). since the lower energy limit 

214 is 200 MeV in our analysis, instead of 100 MeV in the catalog. The relative normalizations of 



|http: //fermi .gsfc .nasa.gov/ssc/data/access/lat/BackgroundModels .html 
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215 the Galactic and isotropic models become 1.02ib0.02 and 1.40ib0.06, respectively, and the fit is 

216 reasonable within the current background model uncertainty. This fit gives a power-law photon 

217 index of Cen A between 200 MeV and 30 GeV of r=2.76 it 0.07 and the flux extrapolated down 

218 to >100 MeV is (2 .06zfc0.20) x 10~^ ph cm~^ s^^ (where errors are statistical only). As noted in 
219 |Abdo_et_a]J (|2009d ) . the spectrum is very steep in comparison to the typical blazars of F = 1.5 — 2 .5 



220 The power-law photon index is consistent with the 3EG result of F = 2.58 it 0.26 (jHartman et al 

221 199S ) . The 3EG flux was reported to be (1.36 ±0.25) x 10 "'^ ph cm s , and have a peak value 

222 of (3.94 lb 1.45) X 10~^ ph cm~'^ s~^ ( Hartman et al.lll999l ). consistent with with the average flux. 



223 We next modeled the region with a radio image of the giant lobe (model B). This analysis is 

224 identical to that described in the lobe paper, and the reader is referre d to it for details. We p resent 

225 a brief description below. We use the WMAP image at 20 GHz from iHardcastle et al.l ()2009l ) . and 

226 eliminate the Cen A core region with a cut radius of 1°. In this analysis, we exclude two point 

227 sources (IFGL J1322.0 — 4515 and IFGL J1333.4 — 4036), which are assumed to be emission from 

228 the lobes. The binned likelihood analysis was performed to extract the flux and spectral indices 

229 for the core and lobes. The relative normalizations of the Galactic and isotropic models become 

230 1.00±0.02 and 1.44ib0.06, respectively. The 7-ray detection in each energy range is significant at a 

231 4(7 level up to the 5.6-10 GeV energy bin for the core region and the spectrum is consistent with 

232 the power-law model. This fit gives a photon index of the core between 200 MeV and 30 GeV of 

233 T=2.67±0.Wstat±0msys and a flux extrapolated down to >100 MeV of {1. 50±0. 25 stat ±0.37 sys) x 

234 1 0^^ ph cm~^ s~^, with statistical and systematic errors reported. Here, we consider the systematic 

235 errors from the effective area, the diffuse model, and WMAP inner cut radius, as described in the 

236 lobe paper. The photon index is almost identical to that of model A, but the flux is somewhat 

237 lower due to some of the core photons from model A being considered as being emitted by the lobes 

238 in model B. The results for model B can be seen in Figure [3l 



3.3. Time Variability 



240 To quantify variability within the ~10 month LAT observation, we generated light curves in 

241 30 and 15 day bins using the unbinned likelihood analysis with gtlike. We performed the analysis 

242 taking into account the lobe emission (i.e.. Model B in § 13. 2p . The power-law normalizations of 

243 the core and background point sources are treated as free parameters, but the photon indices of 

244 all sources and the normalizations of the lobes and the diffuse background models are fixed to the 

245 values obtained in 200 MeV - 30.0 GeV for the whole time region. Figure H^i shows the light curve 

246 of the flux (extrapolated down to > 100 MeV) in 30 day bins. The test results in x^/d.o.f. 

247 = 0.98, and the light curve with 15 day bins gives x'^/d-o.f. = 0.89. These are consistent with 

248 no variability. The time behavior of Cen A is in contra st to large variabil ity of typic al blazars in 



249 the M eV/GeV range, and similar to that of Perseus A (lAbdo et al.ll2009bl ) and M87 (lAbdo et al 
25o l2009dh . 
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Other Contemporaneous Observations 



252 Observations with several different instruments, both on the Earth and in space, were made 

253 during the 10 months of LAT observations presented here. Cen A was observed in the radio as 

254 part of the Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry (TANAMI) 

255 program (Muell er et al. 2009; Oiha et al. 2009). Data were taken with two instr uments on the 



2007 



Kovama et al. 



2007 



256 Swift spacecraft (IGehrels et al.ll2004l') and two ins truments on the Suzaku spacecraft (jMitsuda et al 



Takahashi et al.l 120071 ). A summary of these observations can be found 



258 in Table [U and descriptions are given below. 



4.1. Southern Hemisphere LBA Observations 



260 Cen A was observed with VLBI on 2009 November 27/29, as part of the TANAMI program 

261 using the five antennas of the Australian Long Baseline Array (LBA), the 70m DSS-43 antenna 

262 at NASA's Deep Space Network at Tidbinbilla, Australia, and two trans-oceanic telescopes TIGO 

263 (Chile) and O'Higgins (Antarctica) of the International VLBI Service (IVS) for Geodesy and As- 

264 trometry (the latter two participating at 8.4 GHz, only). The beam size achieved was (0.92 mas x 

265 0.56 mas) at 8.4 GHz and (1.68 mas x 1.25 mas) at 22.3 GHz using natural weighting. These ob- 

266 servations were part of the TANAMI monitoring of a radio and 7-ray selected sample of 65 blazars 

267 at 8.4 GHz and 22.3 GHz with observations approximately every two months. 



268 TANAMI data are correlated on the DiFX software correlator (jDeller et al.l 120071 ) at Curtin 

269 University in Perth, Western Australia. Data inspection and fringe fitting was done with AIPS 

270 (National Radio Astronomy Observatory's Astronomical Image Processi ng System software). The 

271 images were produced by applying the program difmap (ISheDherdlll997l). using the CLEAN algo- 



272 rithm. More details about the data reduction can be found in 



273 Data from the first epoch (November 2009) of TANAMI observations are presented in lOiha et al 



274 (|200g|). Fig.Oincludes the fluxs at 22.3 GHz and 8.4 GHz measured in 2009 November 27/29, respec- 

275 tively. The total flux density, corresponding to the emission distributed over the inner ~ 120 mas 

276 at 8.4 GHz, is 5totai = 3.90 Jy. At 22.3 GHz, a total VLBI flux density of 3.2 Jy is distributed over 

277 the inner ~ 40 mas of the jet, with very little emission on the counterjet side. 

278 Via model fitting, we found a component with an inverted spectrum, which is the brightest 

279 at both frequencies and which we identify with the jet core. The core flux density is 0.92 Jy at 

280 8.4 GHz and 1.54 Jy at 22.3 GHz. The core size is consistently modeled at both frequencies to be 

281 (0.9-1.0) mas x (0.29-0.31) mas at the same position angle of 53-55 degrees (see Ojha et al. 2009). 
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4.2. Suzaku Observations 



2B3 Cen A was observed with Suzaku on 2009 July 20-21, Aug 5-6, and Aug 14-16 with a total 
2B4 exposure of 150 ks, during which time the flux approximately doubled. We utilized data processed 

285 with version 2.4 of the pipeline Suzaku software, and performed the standard data reduction: a 

286 pointing difference of < 1.5°, an elevation angle of > 5° from the earth rim, a geomagnetic cut-off 
2B7 rigidity (COR) of >6 GV. We did not use events from the time the spacecraft entered the South At- 

288 lantic Anomaly (SAA) to 256 s after it left the SAA. Further selection was applied: Earth elevation 

289 angle of > 20° for the X-ray Imaging Spectrometer (XIS), C0R>8 GV and the time elapsed from 

290 the SAA (T_SAA_HXD) of >500 s for the Hard X -ray Detector (HXD ). The XIS response matrices 

291 are created with xisrmfgen and xissimarfgen (jlshisaki et al.l 120071 ). The HXD responses used 

292 here are ae_hxd_pinhxnome5_20070914 . rsp for the PIN and aeJixd_gsohxnom_20060321 .rsp and 

293 aeJixd_gsohxnoni_2007 0424.arf for the Gad olinium Sihcate (GSO) crystal. The "tuned" (LCFIT) 

294 HXD background files (jFukazawa et al.ll2009l ) are utilized. The detailed Suzaku analysis, including 

295 GSO data and time variability, will be reported elsewhere (Y. Fukazawa et al. 2010, in preparation). 

296 The Suzaku data were fit with a single absorbed power-law, which was found to have a spectral 

297 index F = 1.66 it 0.01 with dust absorbing column density Nh = (1.08 ± 0.01) x 10^'^ cm^^. The 



298 flux in the 12 - 76 keV band on 2009 July was (1.23 ± 0.01) x 1 0' 

299 twice the flux measured by Suzaku in 2005 (jMarkowitz et al.l 120071 ). 



ergs ^ cm ^ keV 



about 



4.3. Swift-XBT Observations 

Cen A was obse rved on six days betw een 2009 Jan. 15 - 28 for a total exposure of 22 ksec (see 



301 

302 Tabled]). The XRT (iBurrows et al.ll2005l ) data were processed with the XRTDAS software package 

303 (v. 2.5.1) developed at the ASI Science Data Center (ASDC) and distributed by the NASA High 

304 Energy Astrophysics Archive Research Center (HEASARC) within the HEASoft package (v. 6.6). 

305 Event files were calibrated and cleaned with standard filtering criteria with the xrtpipeline task 

306 using the latest calibration files available in the Swift CALDB. 

307 The XRT dataset was taken entirely in Windowed Timing mode. For the spectral analysis we 

308 selected events in the energy range 2-10 keV with grades 0-2. The source events were extracted 

309 within a box of 40x40 pixels (~94 arcsec), centered on the source position and merged to obtain the 

310 average spectrum of Cen A during the XRT campaign. The background was estimated by selecting 

311 events in a region free of sources. Ancillary response files were generated with the xrtmkarf task 

312 applying corrections for the PSF losses and CCD defects. 

313 The combined January X-ray spectrum is highly absorbed. Hence it was fitted with an ab- 

314 sorbed power-law model with a photon spectral index of 1.98 ±0.05, an intrinsic absorption column 

315 of (9.73 ± 0-26) X 10 ^^ cm~^, in excess of the Galactic value of 8.1 x 10^° cm^^ in that direction 



316 (jKalberla et al.ll2005l ). The average absorbed flux over the 2 — 10 keV energy range is (4.94 it 0.05) x 



317 10 erg cm ^ s ^, which corresponds to an unabsorbed flux of 9.15 x 10 erg cm ^ s 
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318 The XRT spectrum included in the broadband SED was binned to ensure a minimum of 2500 

319 counts per bin and was de-absorbed by forcing the absorption column density to zero in XSPEC, 

320 and applying a correction factor to the original spectrum equal to the ratio of the de-absorbed 

321 spectral model over the absorbed model. 



4.4. Swift-BAT Observations 



323 We used data from the Burst Alert Telescope (BAT) on board the Swift mission to derive a 

324 14-195 keV spectrum of Cen-A contempo rary to the LAT observat ions. The spectrum has been 

325 extracted following the recipes presented in lAjello et al.l (j2008l . l2009bl ) . This spectrum is constructed 

326 by calculating weighted averages of the source spectra extracted over shor t exposures ( e.g. 30 s). 

327 These spectra are accurate to the mCrab level and the reader is referred to lAiello et al.l (j2009al ) for 

328 more details. 



5. SED and Modeling 
5.1. Spectral Energy Distribution 



331 The LAT spectrum of the core of Cen A is shown in Fig. [3l extr apolated into the TeV regime, 

332 along with the HESS spectrum observed between 2004 and 2008 (lAharonian et al.ll2009l ). Also 

333 shown is the HESS spectrum scaled down by its source flux normalization uncertainty. It seems 

334 that the LAT spectrum, with its statistical and systematic errors, extrapolated to higher energies, 

335 is just barely consistent with the HESS spectrum. However, one should keep in mind that the 

336 HESS and LAT spectra presented in this figure are not simultaneous, although the HESS data did 

337 not show any signs of variability. Additionally, 77 absorption makes it unlikely that the HESS and 

338 LAT emission originate from the same region, which is explored below (§ 15. 2p . 

339 Since the cores of many blazars have been shown to be 7-ray loud it is plausible to assume that 

340 the radio core is the source of the central 7 -rays from Cen A . How ever, one should keep in mind that 

341 the error circles of the Fermi and HESS (jAharonian et al.ll2009l ) observations are consistent with 

342 emission from the inner lobes, jet and radio core, so that these other regions could be sources of 

343 7- rays as well. We con s truct the SED for the resolved sub-arcsec and arcsec-scale core as compiled 



344 m 



Meisenheimer et al.l (|2007l ). including their mm/IR/optical observations from 2003-2005. They 
345 have compiled additional points from the 1990s and have applied an ext inction correction of Ay = 9 



346 mag to the optical and I R data. We plot historical data in the X-ray (lEvans et al.ll2004l). hard X 



347 rays 



(iKinzer et al.lll995l : iRothschild et al 



19981 ). and the HESS TeV 7-rays (lAharonian et al 



2006 



Markowitz et al.ll2007l 'l. COMPTEL (ISteinle et al 



20091). The Swift XRT and BAT, as well as 



349 Suzaku data, corrected for Galactic dust as well as dust in NGC 5128, discussed in §[11 were collected 

350 during time intervals which overlap with much of the Fermi-hAT data. Furthermore, we add the 
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351 simultaneous radio data of the TANAMI VLBI jet components. All these are shown in Fig. [5j The 

352 LAT data points in Fig. [5] are from Model B and include statistical errors only. 



5.2. Synchrotron/Synchrotron self-Compton Model 



354 Single-zone synchrotron/synchrotron self-Compton (SSC) models have been v ery successful in 

355 explaining the multiwavelen gth (including 7-ray) emission from BL Lac objects (e.g. jBloom &: Marscher 



19961 : iTavecchio et al.lll998l ). If FRIs are the misaligned counterpart to BL Lacs, one would expect 



357 this model to apply to them as well. In this scenario the low energy, radio through optical emission 

358 originates from nonthermal synchrotron radiation from a relativistically moving spherical homo- 

359 geneous plasma blob, and the X-ray through VHE 7-rays from the Compton scattering of that 

360 synchrotron radiation by electrons in the same blob. The one-zone SSC model has successfully fit 

361 the emissi on from the other Fermi-LAT detected FRIs, Perseus A (NGC 1275: lAbdo et al.ll2009bl ) 

362 and M87 (lAbdo et al.l 2009dl ). and has been successfully applied to previous observations of Cen A 



363 (jChiaberge et al 



2001 



. Here we apply the single-zone SSC model to fit the recent multiwavelength 

364 observations of Cen A, particularly the Fermi-LAT and HESS emission. 

365 One can show (see Appendix [X]) that, on the assumption that all of the emission in the 

366 multiwavelength SED of the Cen A core originates from the same region in a single zone SSC 

367 model, 77 absorption gives the constraint on the Doppler factor 



Sd > 5.3 , 



(1) 



368 where the Doppler factor is 5d = — the bulk Lorentz factor of the jet is Tj = 

369 (1 — I3j)~^^'^, l3jC is the speed of the jet, and 6 = cos~^ /u is the angle of the jet with respect to our 

370 line of sight. Solving for Tj in terms oi 5d, 



371 In order for Tj to be real, the quantity under the radical must be positive, which implies 



5d < 



1 



CSC I 



(2) 



(3) 



372 (e.g., lUrry fc Padovanilll995l ). For Cen A, estimates of 9 vary from 15° to 80° (see section [2]). For 

373 the least constraining value, 9 = 15°, 

5d < 3.8 . (4) 

374 Clearly, the constraints ([1]) and dH are not compatible. Thus, if the radio through Fermi 7-ray 

375 data presented in Fig. [5] are synchrotron and SSC emission originating from the same region of the 

376 jet, then the HESS emission cannot originate from the same part of the jet. Note also that the 

377 HESS emission cannot originate from the same region of the jet, yet be emitted from a different 
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378 mechanism than SSC (say, Compton scattered accretion disk or dust torus radiation) because even 

379 this radiation would be subject to the same 77 attenuation by synchrotron photons. 

380 If the VLBI jet core is assumed to be the origin of the high-energy emission, the TANAMI 
331 core-size measurement can be used to calculate an upper limit on the size of the 7-ray emitting 
382 region of < 0.017 pc = 5.3 x 10^^ cm 3. 1 ). Th is is consistent with the VLBI observations of 



383 iKellermann et al.l (| 19971 ) and iHoriuchi et al.l (j2006l ) , and with a variability timescale of t^, ~ 1 day, 

384 given that the emitting region radius Rb is constrained by the variability time by Rb = 5d c ty. 

385 This variability timescale is consistent with the Suzaku observations, although it is not clear that 

386 the Suzaku X-rays come from the same region as the 7-rays. Using this variab ility timescale and 

387 eqns ()Aip and ()A2p . one gets 5d = 0.6 and B = Q G. More precise modeling (jFinke et al.ll2008l ) 

388 gives the green curve in Fig. [5] with the model parameters in Table [2l This curve demonstrates 

389 the emission can be fit with a Doppler factor of unity. This is consistent with a Lorentz factor 

390 of unity or 7, a degeneracy which can be seen in eqn. ([2|). A stationary, nonrelativistic jet can 

391 explain the entire SEP, ex cept the VHE emission. This fit is similar to the synchrotron/SSC fit by 

392 Meisenheimer et al. (j2007l ) who fit similar data. We further note that a small change in leads 



393 to a large change in the Lorentz factor. This, combined with the uncertainty in the inclination 

394 angle, leads to the fact that the Lorentz factor is not well-constrain ed by modeling. W e also note 

395 that VLBI observations show apparent motion with fij^app ~ 0.1 (jTingav et al.lll998l ). implying 

396 Tj > 1.005, which is also not a particularly strong constraint. 

397 What if the hard X-ray emission originates from thermal Comptonization near the disk, and 

398 not from jet emission? If we assume the rest of the high-energy SED is from the jet, then e^^^ = 1 

399 and /^/'^ = 9 X 10"^^ erg s"^ cm"^, so that eqns (jAT]) and (|A2]) give 5d = 2.4 and 5 = 0.6 G for a 

400 variability timescale of 1 day. More detailed modeling gives the violet curve seen in Fig. [5] with the 

401 parameters in Table [2j The larger Doppler factor needed for this model requires a smaller angle 

402 to the line of sight. The Lorentz factor is again not strongly constrained, and could plausibly be 

403 as high as Vj ~ 8 and still provide a good fit, although this would push the parameters to their 

404 extremes. This model still under-predicts the HESS data. 

405 Jet powers for these models are given in TableEl The proton and pair content of the jet are not 

406 well known, so the total jet power presented in Table [2] is for a pure pair jet, and can be considerd 

407 a lower limit. Even with 10-100 times more energy in ions than leptons, the absolute jet power 

408 is far below the Eddington luminosity for a 10*^ Mq black hole {L^dd = 1-3 x 10'*^ erg s^^). For 

409 the green curve, the parameters assume F = 7. The jet power needed to inflate the giant lobes of 



410 Cen A in their lifetime, as inferred from the radio spectral break, is 10^^ erg s ^ ( Hardcastle et al 



2OO9I ). This value is approximately consistent with the the green curve model presented in Fig. [5l 



412 A possible explanation for the HESS observations is that the TeV emission is produced by 

413 another blob. We show in Fig. [5] (brown curve) that another synchrotron/SSC-emitting blob can 

414 produce the HESS emission without over-producing any of the other multiwavelength data. The 

415 parameters for this blob are in Table [21 although this fit is not unique and many parameter sets 
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416 would fit the HESS data and not contribute at other wavelengths. Other possible origins for the 

417 VHE emission are discussed in § 16.11 



418 5.3. Decelerating Jet Model 

419 Unification models for blazars suggest that FRII galaxies are FSRQs with the jet viewed away 

420 from our line of sight, and similarly FRIs are the parent population of BL Lacs. In this case, 

421 one would expect non-thermal emission from the cores of radio galaxies, de-beamed compared to 

422 blazars. However, the cores of FRIs seem brighter than what is expected from simply de-beamed 

423 emission from BL Lacs, which implies the radio galaxy core emission is from a slower region than 

424 that of BL Lacs, since the beaming angle is related to the bulk Lorentz factor by 9h ~ l/r^- 

425 There are (at least) two possible explanations for this: (1) the jet consists of a faster "spine", 

426 which is responsible for the on-axis blazar emission, inside a slower outer "s heath" , which would b e 



427 responsible for the off-axis emission seen in the cores of radio galaxies (e.g.. IChiaberge et al.ll200d ): 

428 and (2) a decelerating jet model where the on-axis blazar emission is produced by a faster flow 

429 closer to the black hole, and t he off-axis 7-rays seen in radio ga laxies are produced by the slower 



430 flow farther out along the jet (iGeorganopoulos fc: Kazanasll2003l ). 



431 As an example, we provide a flt to the Cen A SED using this decelerating flow, as the blue 

432 curve in Fig. [5j In this model, the high energy emission is due to upstream Compton scattering of 

433 synchrotron photons produced in the slower part of the flow being scattered by energetic electrons 

434 in the faster, upstream part of the flow. The jet starts with a bulk Lorentz factor Tj^max = 5 and 

435 decelerates down to Tj^min = 2 in a length of Z = 3 x 10^^ cm. The injected power law electron 

436 distribution, n{'y) oc 7^^ has an index p = 3.5, and extends from 7mm = 1600, to jmax = 10'', 

437 and the magnetic field at the inlet is = 0.3 G. Jet powers for this model are similar to the 

438 one-zone SSC model fits presented in § 15.21 although this decelerating model fit is particle- rather 

439 than magnetic field-dominated. We also note that the parameters used in this fit are not unique. 



6. Discussion 



6.1. Origin of VHE 7-ray emission 



442 Since the single blob model does not seem to be able to reproduce the broadband SED of 

443 Cen A, could something else be the origin of the VHE 7-rays? We have already shown that 

444 another blob emitting synchrotron and SSC radiation could explain the HESS e mission without 

445 over-producing any of the other data (Fig. [5] brown curve). Lenain et al. (j2008l ) have presented 

446 a model with multiple blobs, moving at different angles to the line of sight from a large opening 

447 angle, to M87 and Cen A (a mong other objects). This model does seem to be able to explain 

448 this SED (jLenain et al.ll2009l ). It has also been suggested that absorbed 7-rays which create e"'"e~ 
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449 pairs, creating an isotropic halo of electrons in the ISM which Compton-sca. tter the host galaxy's 

450 starlight, leading to isotropically-produc ed 7-rays (jStawarz et al.ll2003l . 120061 ) . The HESS data do 

451 seem to match the IStawarz et alj ( 20061 ) predictions Cen A with a galactic magnetic field of 10 

452 /iG. Compton-scattering off of leptons accelerated by the supermassive black hole magnetosphere, 

453 similar to p article acceleration in pulsar s, has been proposed to explain the VHE 7-ray radiation 

454 from M87 ()Neronov &: AharonianI 120071 ) . This could also explain the HESS data from Cen A 

455 separate from the other multiwavelength emission. As we have noted earlier, what we designate 

456 in this paper as the 7-ray "core" actually encon ipasses the radio cor e, jet, and inner lobes of 

457 Cen A. This is also true for the HESS emission. ICroston et al.l (|2009l ) have noted that a shock 

458 front observed in X-rays in the southwest inner lobe could be a source of TeV 7-rays, which seems 

459 consistent with these observations. 



460 Finally, we note that the SED presented here is constructed from non-simultaneous data. 

461 Although Fermi and HESS 7-rays do not show appreciable variability, they could still be vari- 

462 able on longer timescales. Perhaps for a good, simultaneous multiwavelength SED, a one-zone 

463 synchrotron/SSC model could provide a good fit to all of the data. Probably the best way to dis- 

464 criminate between the above models — simple SSC, Compton-scattering emission from a pair halo, 

465 multiple blobs, etc. — is correlated variability between LAT 7-rays and other bandpasses. This 

466 emphasizes the importance of simultaneous multiwavelength data. 



6.2. Origin of UHE Cosmic Rays 



468 The Auger Observatory results indicate some UHECRs could be originating from Cen A (see 

469 § [2]) . The UHECRs could interact with photons at the source and in the extragalactic background 

470 light leading to an observable signature in the HESS band. If the VHE 7-rays originate from cosmic 

471 rays this could account for the discrepancy between HESS and Fermi 7-rays. Based on the green 

472 curve fit presented in Fig. [5] we can analyze whether it is plausible for cosmic rays to originate from 

473 Cen A, keeping in mind that the parameters of that model are not well constrained (§ 3). 

474 The maximum energy to which cosmic rays can be accelerated is limited by the size scale of 

475 the emitting region and the highest energy they can reach before they are cooled. The former 

476 constraint implies that the highest energy a cosmic ray can reach is 



Ez = 4: X 10 
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Z 



B 



6.2 G 



105 S 



6d 



7.0 



eV 



and the latter implies 



Ez = 5.7 X 10^" J- - 




-1/2 



(5) 



(6) 



478 (e.g.. lHillaslll984l : iDermer &: Razzaqudl2010l ). where i;^ ~ 1 is the acceleration efficiency factor, and 

479 e is the elementary charge, Z is the atomic number, and A the atomic mass of the ion. Note 
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480 that these timescales, and all quantities expressed above, are in the frame comoving with the blob, 

481 although for the particular model considered here, 5d=1 so this is not important. 

482 We assume all parameters have values from the green curve model. Thus, it seems for this 

483 model that it is unlikely that protons will be accelerated to energies above ~ 4 x 10^^ eV, although 

484 it is possible for heavier ions to be accelerated this high before they are disintegrated by interacting 

485 with infrared photons from the Cen A core. The threshold energy for photomeson interaction with 

486 peak synchrotron photon s is similar to Ez- Thi s process could create observational signatures from 

487 secondary emission fe.g. . iKachelriefi et al.ll2009l l. as well as convert protons to neutrons, which can 

488 escape as cosmic rays ([Permer et al.l l2009l ) . Again, we note that this result is strongly model- 

489 dependent, and the parameters of this model are not strongly constrained, so this limit should not 

490 be taken too seriously. For example, a small change in the Doppler factor would have little effect 

491 on the model fit, but would require a large change in the bulk Lorentz factor, Vj. A large change 

492 in Tj would significantly affect the highest energy to which particles could be accelerated, as seen 

493 in eqns ([6]) and ([5]). Furthermore, if we are viewing a slower sheath, UHE cosmic rays could be 

494 accelerated in the faster spine beamed away from our line of sight, which could have significantly 

495 different parameters. Acceleration of protons up to 10^'' eV req uires jet powers of Pj > 10'*^ erg s~^, 



496 which may take place in occasional flaring activities in Cen A (jPermer et al. 



2o3). 



497 7. Summary 

498 We have reported on observations of Cen A with the LAT instrument on board the Fermi 

499 Gamma-Ray Space Telescope. This instrument's excellent angular resolution compared to other 

500 7-ray detectors at MeV-GeV energies makes it possible for the first time to separate the lobe and 

501 core emission. The LAT observations have been supplemented with simultaneous observations from 

502 Suzaku, Swift, the Australia Telescope Long Baseline Array, and a variety of non-simultaneous data, 

503 including those from HESS. Our results are as follows: 



1. The LAT-detected core positio n is consistent with C en A's VLBI core (jMa et al.lll998l ) and 
previous EGRET observations (jHartman et al.lll999l ). 



506 2. With 10 months of LAT exposure, we find the core flux > 100 MeV to be (1.50±0.25stat ± 

507 0.37 sys) X 10~^ ph cm~^ s~^ and the spectral ind ex in the 0.2 — 30 GeV range to be r=2.67± 



O.lOstat i O.OSsys, consistent wit h the EGRET (IHartman et al.l Il999l ) and the previously- 
reported 3-month LAT detection ( Abdo et al. 2009a). 



3. 



Extra polated to higher energ ies, the LAT spectrum is barely consistent with the HESS spec- 
trum ([Aharonian et al.ll2009l ) only if the HESS spectrum is lowered in flux by its normalization 



error. 



513 
514 



4. A single zone SSC model can explain all of the multiwavelength emission from the core except 
for the non-simultaneous HESS emission. It is not possible to flt the entire SED, including the 
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515 
516 



HESS emission, with a single zone Compton-scattering model due to internal 77 absorption 
effects. 



5. Modeling results are consistent with suggestions by Chiaberge et al. that we are seeing 7-rays 

from a different origin that we would if were wer e looking down the jet. Th is could be ex- 

plain ed by a spine in sheath (jChiaberge et al. or decelerating jet scenario (jGeorganopoulos &: Kazanad 



20031 ). 
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539 A. 77 Absorption Constraint on the Doppler Factor of Cen A 

540 In the SSC model, the Doppler factor, 6d, and comoving, tangled, isotropic magnetic field 

541 strength, B, may be estimated from the dimensionless peak energy and I'Fi, flux, Cpk and f^f,"" of 

542 the synchrotron and SSC components, respectively, observed in the SED. Assuming the comoving 

543 blob size can be constrained by i?^ = ty6Dc/ (1 + z), this gives 



f^y"- \ /in-10 -1 -2\^^'^ 

Jpk \ / 10 erg s cm \ 



lO^^'^ erg s~^ cm ^ y ^ fpk'^ 



- 20 - 



B = 0.26 G 



1 day 



fSSC 

J pk 



10 erg s ^ cm ^ 




(A2) 



544 (jGhisellini et alJ Il996l ) where is the variabihty timescale and D is the distance to the source. 

545 The Doppler factor, = — (3jfi)]~^ where the Bulk Lorentz factor is Tj = (1 — /3|)~-^/^, (3jC 

546 is the speed of the jet, and 6 = cos~^ /u is the angle of the jet with respect to our line of sight. In 

547 order for 7-rays to escape an emission region, the 77 e~^e~ absorption optical depth, r^^, cannot 

548 be too large. Assuming the uFi, synchrotron flux, fe^^, is given by a broken power law, then for 

549 < 1 for a photon with dimensionless energy e^, this implies 



Sd > 



D 



fsyn 



VIO^; ^1025 cm 
1 day 



(A3) 



t',< 



( Dondi &: Ghisellinilll995l ) where /*^" oc e"^ and A is the index of the synchrotron spectrum below 
the break for 



_i (1 + z)'^ebrk 



< 



26d 



and above the break for 



1 (1 + Z)^eb: 



> 



•rk 



25d 



550 Solving eqn. ()Aip for and inserting this into eqn. 



Sd > 4.4 



syn 
^pk 



2^-^(1 + .)2-2^ (^) 

nsyn 

J pk 



3|), one gets the constraint 
D \ I 



fsyn 



10^^ cm 

V 1/2 



10 1° erg s 1 cm 



(A4) 



1 X 10 



-7 



,ssc 



^pk 



10 erg s 1 cm ^ 



10 erg s ^ cm ^ 
Jssc 

J T) 



I pk 



1/4 



For Cen A, D = 3.7 Mpc = 1.1 x 10^^ cm, and 2; ~ 0. The spectral parameters can be obtained 
from the SED of the core of Cen A (see Fig. [5]) : e^^" 



1.6 X 10" 



,ssc 



pk -0.3,/;r = ^xiU ^"erg 

s~i cm~2, and fp^'^' = 9 x 10"^'^ erg s~i cm"^. Note that here we assume that the X-ray data is 
from the jet; see above. Below the break in the synchrotron spectrum, A ~ 0.5, and above A ~ — 1. 
The highest energy photon bin in the HESS spectrum is e-y = 8 x 10^, so that /*^" = 2 x lO^^'' erg 

s^i cm^2. These values give the constraint 



3x 10- 



-10 



do > 5.3 



551 which is equation ([T]) . 
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Table 1. Summary of multiwavelength observations. 



Instrument 


Observation date 


Exposure time 


Frequency/Energy range 


Australian LBA and IVS 


2009 Nov. 27 


3.6 ks 


22.3 GHz 




2009 Nov. 29 


3.6 ks 


8.4 GHz 


Suzaku XIS 


2009 Jul. 20 - Aug. 16 


150 ks 


0.4 - 10 keV 


Suzaku HXD-PIN 


2009 July 20 - Aug. 16 


150 ks 


10 - 70 keV 


Swift XRT 


2009 Jan. 15 - 28 


22 ks 


0.2 - 10 keV 


Swift BAT 


2008 Aug. - 2009 May 


1.9 Ms 


14 - 200 keV 


Fermi LAT 


2008 Aug. 4 - 2009 May 31 


10 Months 


0.2-30 GeV 



Table 2. Model Parameters. 



Parameter 


Symbol 


Green^ 


Blue2 


Violet^ 


Brown** 


tSuiK l^orentz pactor 




7 n 


5 y 2 


O. 1 


9 n 


Doppler Factor 




1.0 


1.79 — i> 1.08 


3.9 


3.1 


Jet Angle 


e 


30° 


25° 


15° 


15° 


Magnetic Field [G] 


B 


6.2 


0.45 


0.2 


0.02 


Variability Timescale [sec] 


iv 


1.0 X iqs 




1 X 10^ 


1 X 10^ 


Comoving blob size scale [cm] 


Rb 


3.0 X 10^5 


3 X 10^5 


1.1 X lO^'^ 


9.2 X 10*5 


Low-Energy Electron Spectral Index 


Pi 


1.8 


3.2 


1.8 


1.8 


High-Energy Electron Spectral Index 


P2 


4.3 




4.0 


3.5 


Minimum Electron Lorentz Factor 




3 X 102 


1.3 X 103 


8 X 10^ 


8 X 102 


Maximum Electron Lorentz Factor 


'Jmax 


1 X 10^ 


1 X 10^ 


1 X 10* 


1 X 10* 


Break Electron Lorentz Factor 


Ihrk 


8 X 10^ 




2 X 10^ 


4 X lO'^ 


Jet Power in Magnetic Field [erg s"-"^] 




6.5 X 10*3 


1.7 X lO"*! 


2.7 X 10''^ 


4.3 X 10^* 


Jet Power in Electrons [erg s^-"-] 


p 


3.1 X 10"'^ 


3.1 X 10''^ 


2.3 X 10''^ 


7.0 X 10*0 



iSSC Model 

•^Decelerating Jet Model iGeorganopouIos fc KazanajbOOsh 
^SSC Model excluding X-rays 
*SSC Fit to HESS data only 



OO.HESUai 347.8-3751 
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1FGLJ1320.1-4007 
.4-4036 © 1FGLJ1 307.0-4030 

I e 
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• 1FGLJ1 304.3-4352 

1FGLJ1 334.2-4448 * 
©1 FGLJ1 322.0-451 5 
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1FGLJ1 304.0-4622 



1FGLJ1 328.2-4729 



1FGLJ1 305.4-4928 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Fig. 1. — LAT gamma-ray image in the 0.2-30 GeV range in a 14° x 14° region, smoothed by a 
Gaussian with a = 0.3°. The green crosses are the source in the 11 month LAT source hst. Green 
circles are sources considered in the likehhood fitting for model B (see the lobe paper). Red circles 
are additional sources considered in model A. Circle radii represent the semi-major error radius in 
the 11-month catalog. 



Fig. 2. — The LAT localization error c ircles indicated on a 3-color image of Cen A. The image is 
made with the VLA 21 cm image from lCondon et alJ (1996), the optical from Digital Sky Survey 
plates from the UK 48 -inch Schmidt telescope, and an archival Chandra X-ray exposure from 
(jHardcastle et alJ 120071 . OBSID 7797). The 7-ray source is clearly positionally coincident with 
Cen A, enclosing the core, kpc-scale jet, and most of the radio lobes. 
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Fig. 3. — Spectrum of the Cen A core from differential fluxes derived for successive energy ranges 
from model B (black circles). The black bowtie indicates the best fit 0.1 - 30 GeV LAT flux and 
r with statistical errors only, while the green bowtie indicates this with systematic errors as well. 
The LAT spectrum is extrap olated into the HESS energy range (dashed lines). The HESS data 
from lAharonian et al.l (|2009l ) are shown (red squares) and the HESS data shifted to lower flux by 
their statistical and systematic normalization error (blue squares). The latter are also shifted in 
energy by 10% for clarity. 
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Fig. 4. — (a) Fermi-LAT light curve of Cen A without considering lobe emission (Model A) in 30 
day bins, with (b) simultaneous lightcurve from Swift-Ji AT (14 day bins). 
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Fig. 5. — The SED of the Cen A core with model fits. Colored symbols are observations between 
August and May 2009, the epoch of the LAT observations. These include observations of, from 
low to high frequency: the TAN AMI VLBI (red squares), Swift-XI{T (red crosses), Suzaku (brown 
circles), Swift-B AT ( red circles), and Fermi-L AT (re d diamonds). Black sy mbols are archival data, 
(jMarconi et al.ll2000l ) including HESS observations (jAharonian et al.l 12003 ). Curves are model fits 
to nuclear region of Cen A. The green curve is a synchrotron/SSC fit to the entire data set. The 
dashed green curve shows this model without 77 attenuation. The violet curve is a similar fit but 
is designed to under fit the X-ray data, and the brown curve is designed to fit the HESS data while 
not over-producing the other data in the SED. The blue curve is the decelerating jet model fit 
(jGeorganopoulos fc Kazanasll2003l ). See Table [2] for the parameters of these model curves. 



